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Abstract 

We study the effect of transport processes (diffusion and free-streaming) on a 
collapsing spherically symmetric distribution of matter in a self-similar space- 
time. A very simple solution shows interesting features when it is matched 
with the Vaidya exterior solution. In the mixed case (diffusion and free- 
streaming), we find a barotropic equation of state in the stationary regime. 
In the diffusion approximation the gravitational potential at the surface is 
always constant; if we perturb the stationary state, the system is very sta- 
ble, recovering the barotropic equation of state as time progresses. In the 
free-streaming case the self-similar evolution is stationary but with a non- 
barotropic equation of state. 
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I. INTRODUCTION 



The formation of compact objects is usually preceded by an epoch of radiative collapse 
[1]. One of the main difficulties in the study of these systems is that there is no rehable 
information about the equation of state in the central regions of superdense stars, such as 
neutron stars, and this leads to assumptions of a very general natTirc [2]. 

Few exact solutions to the Einstein equations are relevant to gravitational collapse. For 
this reason, new collapse solutions are very useful, even if they are simplified ones [3]. It 
is well known that the field equations admit homothetic motion [3-6]. Applications of 
homothetic similarity range from modeling black holes to producing counterexamples to the 
cosmic censorship conjecture [7-14]. 

The natural formulation of the collapse problem is that of initial values, with an a priori 
defined equation of state. However, there are only a few solutions to the Einstein field 
equations with a well defined equation of state. It has been shown that the only perfect 
fluid equation of state compatible with self-similarity is the barotropic one [4] . 

Emission of photons or neutrinos is a typical process in the evolution of massive stars. 
The only plausible mechanism to deliver almost all the binding gravitational energy, during 
the collapse toward a neutron star, is that of neutrino emission [15]. It seems clear that the 
free-streaming process is associated with the initial stages of the collapse, while the diffusion 
approximation becomes valid toward the final stages. 

In this paper, we explore self-similar gravitational collapse. For such an assumption, 
we recast the geometrical variables so as to have an explicit radial dependence. The fiuid 
is considered with heat fiow (the diffusion approximation) or free propagation in the radial 
direction (free-streaming) . The interior solution is matched to an exterior ( Vaidya) solution 
by means of the Darmois-Lichnerowicz conditions. 

We discuss neither the micropro cesses that produce free-streaming, nor the tempera- 
ture distribution during diffusion. For the latter another approach is necessary to avoid 
pathological behavior [16]. 



For the geometric description of the interior of the distribution we use the radiation 
metric in the spherical Bondi form [17] 



where (3 and V are functions of u and r. Here m is a time-like coordinate, r is a null 
coordinate {grr = 0), that is, r > is an affine parameter along the null generators of 
u — constant null hypersurfaces, and 9, (p are the usual angular coordinates; we are using 
geometrized units {c — G— 1). 

The hydrodynamic scenario, as viewed by a local Minkowskian observer comoving with 
the fiuid (with velocity —u), consists of an isotropic fiuid of density p and pressure p, 
unpolarized energy density e and heat fiux q traveling in the radial direction. Therefore, for 



II. THE FIELD EQUATIONS AND MATCHING 
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this comoving observer, the covariant energy-momentum tensor is 
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Note that the velocity of matter in the Bondi coordinates is 

dr V u 
du r 1 — uj 

We can write the Einstein field equations as [18] 



(2) 



[p + puj"^ + 2ujq + e{l+ujf]=m,r- T-^ (3) 



[{p-ujp)-{l-uj)q]=m^r, (4) 



(l + o;) 

27rr(l -a;)(p + p-2g) 
(l-2m/r)(l + u;) 



= (5) 



Sttp = -2e-2^/3,„, + [3/3,,(l - 2m,,) - m,,,]/r + (1 - 2m/r)(2/3,,, + 4/3^ - /3_^/r), (6) 

where the comma subscript represents partial difi'erentiation with respect to the indicated 
coordinate and m is the Bondi mass defined by m = (r — F exp(— 2/3))/2. 

The exterior space-time is described by the Vaidya radiating metric [19]. In order to 
match this to the interior solution, we use the Darmois-Lichnerowicz conditions. These are 
equivalent to the continuity of the functions {3 and in across the boundary of the sphere, 
and to the condition [20,21] 

[-/3,„e2^ + (1 - 2rh/r)(5,r - m,,/(2r)]„ = , (7) 

which is equivalent to Pa = qa- The subscript a indicates that the quantity is being evaluated 
at the surface r = a{u). 

III. SPHERICAL AND SELF-SIMILAR INTERIOR SOLUTIONS 

Self-similarity is invariably defined by the existence of a homothetic Killing vector field 
[4]. A homothetic vector field on the manifold is one that satisfies £^g =2ng on a local 
chart, where n is a constant on the manifold and £ denotes the Lie derivative operator. If 
n 7^ we have a proper homothetic vector field and it can always be scaled to have n — 1; 
if n = then ^ is a Killing vector on the manifold [22-24]. So, for a constant rescaling, ^ 
satisfies —2g and has the form ^ = A(m, r)du + X{u, r)dr. If the matter field is a perfect 
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fluid, the only equation of state consistent with £^g =2g is a barotropic one [4] . The homo- 
thetic equations reduce to ^{X) = 0, ^{Y) = 0, A = r and A = A(m), where X = fh/r and 
Y = Ae^^/r. Therefore, X = X{() and Y = Y{Q are solutions if the self-similar variable is 
defined as ( = r exp(— J du/A). Here we assume that X — CiC,^ and that Y — C2C^ where 
Ci, C2, k and I are constants. This power-law dependence on C, is based on the fact that 
any function of C, is solution of =2g. As we shall see, this simplifying assumption is not 
devoid of physical meaning. Demanding continuity of the first fundamental form we get the 
following metric solutions: 

m = rh,{T/af-^\ (8) 



e2^ = (r/ay+^ (9) 

Condition (7) then implies that the local radial velocity uj is determined at the surface in 
terms of the gravitational potential fha/a and the parameters k and I: 

..= i-^i^^. (10) 

ma/ a (1 + k) 

Equations (2) and (3) evaluated at the surface constitute the system of (ordinary differential) 
equations in a{'u) and rha{u), to be integrated while taking into account Eq. (10). Thus, the 
dynamics at the surface is completely determined if we establish how energy is exchanged 
with the exterior. In fact, we have found that self-similarity determines the luminosity 
profiles [25,26]. 



IV. TRANSPORT PROCESSES 

In order to explore the effect of the transport processes on self-similar gravitational 
collapse, we consider below a combination of the diffusion and free-streaming mechanisms, 
and the action of each one separately. In particular we discuss whether the barotropic 
equation of state holds in each case. 



A. Mixed 

If the transport mechanism is mixed, we require additional information. Only in this 
case do we suppose orthogonality between the four-velocity and the homothetic vector. This 
condition has been employed to obtain static solutions [27] and it also establishes a rela- 
tionship between the homothetic vector and the equation of state [28]. Another meaning 
attributable to the orthogonality condition is that the group, generated by the homothetic 
vector, is acting upon the three-space comoving with the observer. The referred condi- 
tion, together with £^g =2g, can thus be seen as a covariant definition of self-similarity 
[29]. Therefore, from the orthogonality condition we obtain the temporal component of the 
homothetic vector: 

A = m„(l + 0/(l + A;). (11) 
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The radiation flux at the surface is then determined from the symmetry equations. The 
heat flux at the surface is also determined from equation (6) or, equivalently, from T|!^ = 
evaluated at the surface. Now, feeding back (10) and (11) into the symmetry equations, 
we find algebraically by means of symbolic manipulation with REDUCE, two restrictions 
which allow us to satisfy the symmetry equations: (i) k — —I; (ii) k — [{1 — 2ma/a){l + 2) — 
l]/(2ma/a). Restriction (ii) leads to uja = —1 (through Eq. (10)), that is, the fluid collapses 
at light speed; we must reject this possibility. Under restriction (i), k = —l, the homothetic 
symmetry is preserved at all points of the space-time, when we integrate numerically (using 
the Runge-Kutta method) the system of equations at the surface. Specifically, for A; = 
the physical variables are well behaved and the equation of state is barotropic in the whole 
sphere; that is, the ratio p/p is constant at all points of the material (see Figure 1). All 
shells collapse with the emission of energy. The pressure, density, heat flow and radiation 
flux are stationary inside the sphere. The matter velocity decreases toward the center and 
is constant at the surface. Observe that p ^ (3/5)p, as we expect, because the heat flow 
diminish the effective gravitation or, in others words, the equation of state is softened by 
the diffusive process [30,18]. For k ^ the equation of state is not longer barotropic and 
the physical variables are not stationary, although the interior space-time is self-similar. 



B. Diffusion 



In the diffusion approximation (e = 0) we again obtain the heat flux at the surface 
from equation (6) or T^.^ — evaluated at the surface. From the symmetry equations 
evaluated at the surface, we deduce that the gravitational potential at the surface, rha/a, is 
a function of only k and /. The surface equations thus reduce to one differential equation 
for a (or rha). This situation leads us to four possible restrictions: (i) k = —I; (ii) / = — 1 
which is equivalent to a;,! = 1 (see Eq. (10)), that is, the fluid explodes at light speed; 
(iii) a complicated polynomial of degree three in k with coefficients depending on /; (iv) a 
complicated polynomial of degree seven in k with coefficients depending on I. We solved 
analitycally the polynomial of degree three (using REDUCE), obtaining one real and two 
complex solutions. All these must be rejected because they do not have physical meaning 
for a wide interval of numerical values of /. After solving numerically the polynomial of 
degree seven, we found that the physically acceptable models are similar to those that 
emerge from restriction (i), k — —I. In any case, the gravitational potential at the surface 
is constant {rha/a depending only on k) although the distribution collapses. If we perturb 
this stationary state, we find numerically that the system recovers the barotropic equation 
of state only for k = (see Figure 2). Perturbations consist of enhancing or diminishing the 
gravitational potential at the surface, while preserving the values of k. Therefore, we must 
integrate numerically (using the Runge-Kutta method again) the two perturbed differential 
equations at the surface. It is important to stress that the symmetry equations are satisfied 
everywhere and at all times when we perturb the system. 
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C. Free— streaming 



In this transport process {q = 0) the boundary condition reduces to = 0. We obtain 
again a gravitational potential which is constant at the surface (from T{!^ evaluated at the 
surface). Therefore, the luminosity profile and the temporal component of the homothetic 
vector are determined from one of the symmetry equations. No restrictions appear on the 
parameters k and I, and no barotropic behavior is found, at least for the cases examined. 

V. CONCLUSION 

We have studied the effect of diffusion and free-streaming on self-similar gravitational 
collapse using the solutions to the Einstein equations given by (8) and (9). We find that the 
equation of state is barotropic only when diffusion occurs and the parameter k is zero. In 
general {k ^ 0) the equation of state is not barotropic, although self-similarity holds every- 
where. Thus for non-pcrfcct fluids self-similarity is connected with a barotropic equation of 
state only under special conditions. The heat flow seems to be crucial for this result. Vis- 
cosity in self-similar distributions has also been considered previously and the same result 
emerges; that is, the equation of state can be barotropic [26]. Wc have found relativistic 
and self-similar examples compatible with a nonbarotropic equation of state [31], at least 
within our chosen simple solutions and the orthogonality condition between the generator 
of the group and the four-velocity. It would be of interest to attempt to relax some of our 
suppositions in order to investigate how general our conclusions are. 
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FIG. 1. Evolution of the physical variables for the mixed transport process with A; = 0. In each 
graph the curves correspond (from the uppermost to the lowermost) to ratios: r/a(0) = 0.2, 0.4, 
0.6, 0.8 and 1.0, respectively. The ratio pj p is the same at any point of the material. 
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FIG. 2. The ratio p/p as a function of retarded time in the diffusion approximation with k 
considering perturbations of the stationary regime. 



